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Abstract 

A search is presented for hidden-sector bosons, y, produced in the decay 
> A*(892)°X, with A* (892)° —)• and y —)• The search is 

performed using pp-collision data corresponding to 3.0 fb“^ collected with the 
LHCb detector. No significant signal is observed in the accessible mass range 
214 < m{x) < 4350 MeV, and upper limits are placed on the branching fraction 
product )• A*(892)°x) x R(x—)• as a function of the mass and lifetime 

of the X boson. These limits are of the order of 10“° for x lifetimes less than 100 ps 
over most of the m(x) range, and place the most stringent constraints to date on 
many theories that predict the existence of additional low-mass bosons. 
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Interest has been rekindled in hidden-sector theories (^, motivated by the current lack 
of evidence for a dark matter particle candidate and by various cosmic-ray anomalies [2]-[^ . 
These theories postulate that dark matter particles interact feebly with all known particles, 
which is why they have escaped detection. Such interactions can be generated in theories 
where hidden-sector particles are singlet states under the Standard Model (SM) gauge 
interactions. Coupling between the SM and hidden-sector particles may then arise via 
mixing between the hidden-sector held and any SM held with an associated particle that is 
not charged under the electromagnetic or strong interaction (the Higgs and Z bosons, the 
photon, and the neutrinos). This mixing could provide a so-called portal through which a 
hidden-sector particle, y, may be produced if kinematically allowed. 

Many theories predict that TeV-scale dark matter particles interact via GeV-scale 
bosons (9}jII] (c = 1 throughout this Letter). Previous searches for such GeV-scale 
particles have been performed using large data samples from many types of experiments 
(see Ref. for a summary). These searches have placed stringent constraints on the 
properties of the hidden-sector photon and neutrino portals; however, the constraints on 
the axial-vector and scalar portals are significantly weaker. 

One class of models involving the scalar portal hypothesizes that such a y held was 
responsible for an inhationary period in the early universe 13 , and may have generated 


the baryon asymmetry observed today The associated inhaton particle is expected 

to have a mass in the range 270 < U7.(y) ;< 1800 MeV 13 . Another class of models 


invokes the axial-vector portal in theories of dark matter that seek to address the cosmic- 
ray anomalies, and to explain the suppression of charge-parity (CP) violation in strong 
interactions |^. These theories postulate an additional fundamental symmetry, the 
spontaneous breaking of which results in a particle called the axion |17 . To couple the 


axion portal to a hidden sector containing a TeV-scale dark matter particle, while also 
explaining the suppression of CP violation in strong interactions. Ref. 
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proposes an 

axion with 360 < m(y) < 800MeV and an energy scale, /(y), at which the symmetry 
is broken in the range 1 < /(y) < 3TeV. A broader range of m(y) and /(y) values is 
allowed in other dark matter scenarios involving axion(-like) states [l9| - [^ . 

This Letter reports a search for a hidden-sector boson produced in the decay P° —)■ K*^Xi 
with y —)■ and K*^ —)■ K~^n~ (throughout this Letter, K*^ = A*(892)° and the 

inclusion of charge-conjugate processes is implied). Enhanced sensitivity to hidden-sector 
bosons arises because the 6 —)■ s transition is mediated by a top quark loop at leading 
order (see Fig. [^. Therefore, a y boson with 2m(/i) < m(y) < m(P°) — m(A*°) and 
a sizable top quark coupling, e.g. obtained via mixing with the Higgs sector, could be 
produced at a substantial rate in such decays. The P° —?• A*°y decay is chosen instead 
of —)■ A+y, since better y decay time resolution is obtained due to the presence 

of the A’''7r“ vertex, and because there is less background contamination. The data 
used correspond to integrated luminosities of 1.0 and 2.0 fb~^ collected at center-of-mass 
energies of ^/s = 7 and 8TeV in pp collisions with the LHCb detector. This is the first 
dedicated search over a large mass range for a hidden-sector boson in a decay mediated 
by a 6 —)■ s transition at leading order, and the most sensitive search to date over the 
entire accessible mass range. Previous limits set on y boson production in such decays 
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Figure 1: Feynman diagram for the decay )• K*^x-: with y—)> . 


have either focused on a limited mass range |2^, or have been obtained from more general 
searches for long-lived particles 


23 


The LHCb detector is a single-arm forward spectrometer covering the pseudorapidity 


range 2 < r/ < 5, designed for the stndy of particles containing 6 or c qnarks |24,25 


The detector includes a high-precision charged-particle tracking system for measuring mo¬ 
menta |26p7l ; two ring-imaging Cherenkov detectors for distinguishing charged hadrons 


a calorimeter system for identifying photons, electrons, and hadrons; and a system for 
identifying mnons 29 . The trigger consists of a hardware stage, based on information 


from the calorimeter and mnon systems, followed by a software stage, which applies a 
fnll event reconstrnction 30 . The selection of —)• K*^x candidates in the software 
trigger reqnires the presence of a vertex identified by a mnltivariate algorithm as 
being consistent with the decay of a fe hadron. Alternatively, candidates may be selected 
based on the presence of a displaced dimuon vertex, or the presence of a muon with large 
transverse momentnm (pt) and large impact parameter (IP), dehned as the minimnm 
track distance with respect to any pp-interaction vertex (PV). Only tracks with segments 
reconstructed in the hrst charged-particle detector, which surrounds the interaction region 
and is about Im in length j^, can satisfy these trigger reqnirements; therefore, the y 
boson is reqnired to decay well within this detector. In the simnlation, pp collisions are 
generated following Refs. 32 -35 , and the interactions of the outgoing particles with the 


detector are modelled as in Refs. 36,37 


A search is conducted, following Ref. 38 , by scanning the distribution for an 

excess of y signal candidates over the expected background. In order to avoid experimenter 
bias, all aspects of the search are fixed withont examining those —)■ K*^x candidates 
which have an invariant mass consistent with the known B^ 


mass 
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The step sizes 

in m(y) are where fi~)] is the dimnon mass resolntion. Signal 

candidates satisfy — m{x)\ < 2cr[m(p+/i“)], while the background is estimated 

by interpolating the yields in the sidebands starting at 3(j[m(/i+/r“)] from m(y). With 
m{K~^7r~^~) constrained 40| to the known B^ mass, is less than 8MeV 


over the entire ) range, and is as small as 2 MeV below 220 MeV. The statistical 

test at each m(y) is based on the prohle likelihood ratio of Poisson-process hypotheses with 
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and without a signal contribution 41 


The uncertainty on the background interpolation 
(see Ref. 


38 for details). 


is modeled by a Gaussian term in the likelihood 

The /i'*'/!” decay vertex is permitted, but not required, to be displaced from the 
—)■ K*^x decay vertex. Two regions of reconstructed dimuon lifetime, 


are 


defined for each m{x) considered in the search: a prompt region, |r(/i+/i“)| < 3cT[r(/i’''/i“)], 
and a displaced region, ^~)]. The lifetime resolution is about 0.2ps 

for > 250MeV, and Ips near 2m(/i). The joint likelihood is formed from the 

product of the likelihoods for candidates populating the prompt and displaced regions, 
since no assumption is made about t(x). Narrow resonances are vetoed by excluding the 
regions near the u, 0, J/'ip, '0(25') and ■0(3770) resonances. These regions are removed in 
both the prompt and displaced samples to avoid contamination from unassociated dimuon 
and K*^ resonances. 

The branching fraction product —)■ 77*°x(/r+/r“)) = B{B^ —)■ K*^x) x 

is measured relative to B{B^ ^ where the normalization sample is taken 

from the prompt region and restricted to 1.1 < m?{fB'iJL~) < 6.0 GeV^. This normalization 
decay is chosen since the detector response is similar to that for the B^ —)■ K*^x decay, 
and because the hidden-sector theory parameters can be obtained from the ratio B{B^ ^ 


K*^x{l^~^l^~))/ 13{B^ —)■ with reduced theoretical uncertainty. Gorrelations 

between the yields of a possible signal in the prompt 1.1 < fj.~) < 6.0 GeV^ region 

and the normalization decay are at most a few percent and are ignored. 

with the exception that the K*^ and dimuon 


The selection is similar to that of Ref. 42 


candidates are not required to share a common vertex. Signal candidates are required to 
satisfy a set of loose requirements: the B^, K*^ and x decay vertices must all be separated 
from any PV and be of good quality; the B^ IP must be small, while the IP of the kaon, 
pion and muons must be large; the angle between the B^ momentum vector and the vector 
between the associated PV and the B^ decay vertex must be small; and the kaon, pion 
and muons must each satisfy loose particle identihcation requirements. Candidates are 
retained if m{K~^7r~) is within 100 MeV of the known K*^ mass 
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A multivariate selection is applied to reduce the background further. The uBoost 
algorithm is employed to ensure that the performance is nearly independent of m{x) 
and r(x)- The inputs to the algorithm include pt{B^), various topological features of the 
decay, the muon identihcation quality, and an isolation criterion designed to suppress 
backgrounds from partially reconstructed decays. Data from the high-mass sideband, 
150 < m{K~^7T~fi~) — m{B^) < 500 MeV, are used to represent the background in the 
training, while simulated samples generated with m[x) values of 214, 1000, and 4000 MeV, 
and r(x) large enough to populate the full reconstructible region, are used for the signal. 
The multivariate selection requirement is determined by maximizing the hgure of merit of 
Ref. 45 for hnding a signal with a signihcance of hve standard deviations. This results in 


a signal selection efficiency of 85% with a background rejection of 92% on average. The 
uBoost algorithm is validated using ten additional signal samples generated with various 
other m{x) and t{x) values. The performance is consistent for all samples. 

Peaking backgrounds that survive the multivariate selection are vetoed explicitly. A 
small number of )■ decays are removed by rejecting K^7i~ candidates 
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m(K*7t jX+H ) [MeV] 

Figure 2: Invariant mass spectrum with fit overlaid for all prompt ^ candidates 

with 1.1 < < 6.0GeV^. 


that are consistent with the decay 0—)- K^K~ if the tt~ is assnmed to be a misidentihed 
K~. A similar veto is applied that removes about 250 —)■ decays. Candidates 

are also rejected if the dimuon system is consistent with any of the following decays: 

7r~^n~, where the pions decay in flight to muons; ^ where the pion decays 
in flight and the proton is misidentihed as a muon; and )■ where the kaon and 

pion decay in hight. All other particle-misidentihcation backgrounds are negligible. 

Figure shows the K^t:~ jjA ijl~ mass distribution for all prompt candidates that satisfy 
the full selection in the region 1.1 < < 6.0 GeV^. An unbinned extended maxi¬ 
mum likelihood ht is performed to obtain the —>■ pT yield. The signal model is 

obtained from data using the subset of prompt candidates with m{pApi~) in the J/ip region, 
where the background is 0(10“^). A small correction, obtained from simulation, is applied 
to account for the difference in signal shape expected in the 1.1 < m?{pApi ) < 6.0 GeV^ 
region. The background model is an exponential function. Several alternative background 
models are considered, with the largest shift observed in the signal yield (1%) assigned as 
a systematic uncertainty. The S'-wave fraction (he. not a K*^ meson) of the Ktt system 
within the selected iFvr mass range is (4 ± 4)% [^. The yield of the normalization mode 
is N{B^ — K*^pAp~) = 506 ± 33, where the uncertainty includes both statistical and 
systematic contributions. 

Probability density functions, obtained from the data using splines, are used to generate 
simulated data sets under the no-signal hypothesis from which the global signihcance 
of any y signal is obtained j^. For this the data are collected in the prompt region 
into wide bins with a width of 200 MeV, and into a total of three bins in the displaced 
region. Simulated events show that the presence of a narrow y signal anywhere in the 
m(y)-r(y) plane, whose local significance is 5 ct, would not produce a significant excess in 
these wide-binned data. 

Figureshows the m{pApr) distributions in both the prompt and displaced regions 
for candidates whose invariant mass is within 50 MeV of the known B^ mass. The most 
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Figure 3: Distribution of in the (black) prompt and (red) displaced regions. The shaded 

bands denote regions where no search is performed due to (possible) resonance contributions. 
The tpi^S) and '0(3770) peaks are suppressed to better display the search region. 


significant local excess occnrs for m{x) = 253 MeV, where in the prompt region 11 (6.2) 
candidates are observed (expected), while the displaced region contains a single candidate 
which is the only displaced candidate below m{u:). The p-value of the no-signal hypothesis 
is about 80%, showing that no evidence is found for a hidden-sector boson. 

To set upper limits on i3(5°—)■ various sources of systematic uncertainty 

are considered. The limits are set using the prohle likelihood technique 46 , in which 


systematic uncertainties are handled by including additional Gaussian terms in the likeli¬ 
hood 38 . Since no contamination from the u or (j) resonance is found in the displaced 


region, upper limits are set in these m{x) regions for r(x) > Ips. 

Many uncertainties cancel to a good approximation because the signal and normalization 
decays share the same hnal state. The dominant uncertainty on the efficiency ratio 
e(i?°—)■ if'*°x(/i’''/r“))/e(i?°—)■ which is taken from simulation, arises due to its 

dependence on The simulation is validated by comparing r(7r+7r“) distributions 

between )• J/'ip decays reconstructed in simulated and experimental data in 
bins of momentum. The distributions in data and simulation are consistent in each 
bin, and the per-bin statistical precision (5%) is assigned as systematic uncertainty. 

The uncertainty on the efficiency for a signal candidate to be reconstructed within 
a given m(/r+p“) signal window, due to mismodeling of is determined to 

be 1% based on a comparison of the 0/0 peak between )■ J/0decays in 
simulated and experimental data. A similar comparison for (j[r(/r+/i“)] shows that the 
uncertainty on the fraction of signal candidates expected to be reconstructed in the prompt 
and displaced regions is negligible. Finally, the efficiency for the normalization mode 
is determined using the measured angular distribution 47 , which is varied within the 


uncertainties yielding an uncertainty in the normalization-mode efficiency of 1%. The 
individual contributions are summed in quadrature giving a total systematic uncertainty 
of 
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Figure 4: Upper limits at 95% CL for (left axis) B{B^^ ^ iL*°iU+/i“), with 

5°—)■ in 1.1 < < 6.0 GeV^, and (right axis) B{B^ ^ K*^x{b~^ l^~))- The 

sparseness of the data leads to rapid fluctuations in the limits. Excluding the region near 2m(//), 
the relative limits for r < 10 ps are between 0.005-0.05 and all relative limits for r < 1000 ps are 
less than one. 


The spin of the hidden-sector boson determines the angular distribution of the decay 
and, therefore, affects the efficiency. The upper limits are set assuming spin zero. For 
a spin-one y boson produced unpolarized in the decay, the sensitivity is about 10-20% 
better than for the spin-zero case. The dependence on the polarization in the spin-one 
case is provided as supplemental material to this Letter [48] . 

Figure|^shows the upper limits on —)■ (/r+/r“)), relative to —)■ iF*°yU,+/i“) 

in the 1.1 < m^(p.+/i“) < 6.0 GeV^ region, set at the 95% confidence level (CL) for several 
values of T(y); limits as functions of r(y) are provided as supplemental material to this 
Letter. The limits become less stringent for r(y) > 10 ps, as the probability of the 
y boson decaying within the hrst charged-particle detector decreases. The branching 
fraction B{B^ ^ iF*°/i+/i“) = (1.6 ± 0.3) x 10“^ is used to obtain upper limits on 


B{B^ —^ B~))i which are also shown in Fig. Due to the uncertainty on the 

normalization-mode branching fraction, there is not a one-to-one mapping between the 
two axes in the hgure; however, the absolute limits shown are accurate to about 
Figure 1^ shows exclusion regions for the DFSZ 49,50 axion model of Ref. 


20 set 


in the limit of large ratio of Higgs-doublet vacuum expectation values, tan /5 > 3, for 


charged-Higgs masses m{h) = 1 and lOTeV (this choice of restricted parameter space 
is made for ease of graphical presentation). The constraints scale as log {m{h)/ TeV) for 
m{h) > 800 GeV. The branching fraction of the axion into hadrons varies greatly in 
different models. Figure]^ shows the results for two extreme cases: f3(y—)■ hadrons) = 0 
and 0.99. While f3(y—?• /i^/i“) is 100 times larger when i3(y—)■ hadrons) = 0, r(y) is also 
larger, which results in the model probing the region where the upper limits are weaker. 
The constraints are loose for m(y) > 2m(r), since the axion preferentially decays into 
T^T~ if kinematically allowed; otherwise the exclusions reach the PeV scale. 
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Figure 5: Exclusion regions at 95% CL: (left) constraints on the axion model of Ref. [M; (right) 
constraints on the inflaton model of Ref. . The regions excluded by the theory |51] and by 


the CHARM experiment 52 are also shown. 


Figure]^ also shows exclusion regions for the inflaton model of Ref. j^, which only 
considers m{x) < 1 GeV. The branching fraction into hadrons is taken directly from 
Ref. and, as in the axion model, is highly uncertain but this does not greatly affect 
the sensitivity of this search. Constraints are placed on the mixing angle between the 
Higgs and inflaton helds, 9, which exclnde most of the previonsly allowed region. 

In summary, no evidence for a signal is observed, and upper limits are placed on 
i3(i?°—)■ K*^x) X This is the first dedicated search over a large mass range 

for a hidden-sector boson in a decay mediated by a 6—)■ s transition at leading order, and 
the most sensitive search to date over the entire accessible mass range. Stringent constraints 
are placed on theories that predict the existence of additional scalar or axial-vector helds. 
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Supplemental Material 

The limits reported in the Letter assume a spin-zero hidden-sector boson. To convert these 
into limits for a spin-one boson, the ratio of efficiencies for the spin-one to spin-zero cases 
must be accounted for. Determining this ratio involves integrals of the form 

f /j(f1)e(D, m^(/i+/r“))df1 
/ fic(Q)€(Q,m^ju+ju-))dQ’ 

where D = (see Appendix A of Ref. [40] in the Letter for details on the 

angular basis), is the efficiency, /j(D) are functions of the angles, and 

/ic(D) = cos^ 9 k. Figure shows the values for 

= sin^ 9 k, 

/2s(D) = sin^^if cos26'£, 

/2c(^) = cos^ 9k cos 29£. 

All other integrals, each of which has a value of zero in the absence of inefficiency, have 
values 0(0.01). Therefore, the following terms in the general angular distribution can be 
ignored when determining the limits: 


fsm 

= sin^ 9k sin^ 9i cos 2(/, 

um 

= sin 29k sin 20^ cos 0, 

fm 

= sin 29k sin 9^ cos 0, 


= sin^ 9k cos 9i, 


= cos^ 9 k cos 9^, 

/7(fi) 

= sin 29k sin 9^ sin </, 

fsm 

= sin 29k sin 29^ sin 0, 

fdn) 

= sin^ 9k sin^ 9i sin 20. 


Figure shows an example efficiency ratio for the case of a spin-one boson produced 
unpolarized in the decay. Since the j = 3,4,... 9 terms integrate to approximately zero, 
this same curve applies for any theory that predicts that the longitudinal polarization 
fraction of the K*^ is = 1/3. 
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m(ju-^/u ) [GeV] 



Figure 6: (left) Integral values for Is, 2s and 2c relative to the value for Ic (see text for details). 
The dashed lines show the values in the absence of inefficiency, (right) Ratio of the efficiency for 
an unpolarized spin-one boson to that of a spin-zero boson. 



m(jj.+/j. ) [MeV] 

Figure 7: Upper limits at 95% CL for (left axis) R(i?° —)• ^ 

with )• in 1.1 < < 6.0 GeV^, and (right axis) ^ 

Same as Fig. 4 in the Letter but including the r = 0 and 1 ps limits. 
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